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Abstract In this paper, we modify the Newton-Secant method with third 
order of convergence for finding multiple roots of nonlinear equations. Per 
iteration this method requires two evaluations of the function and one eval¬ 
uation of its first derivative. This method has the efficiency index equal to 
33 « 1.44225. We describe the analysis of the proposed method along with 
numerical experiments including comparison with existing methods. Moreover, 
the dynamics of the proposed method are shown with some comparisons to 
the other existing methods. 

Keywords Multi-point iterative methods • Newton-Secant method • multiple 
roots • Basin of attraction. 


1 Introduction 

Solving nonlinear equations based on iterative methods is a basic and ex¬ 
tremely valuable tool in all fields in the science as in economics, engineering 
and physics. The important appearance related to these methods are order 
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of convergence and number of function evaluations. Therefore, it is favorable 
to attain as high as possible convergence order with fixed number of func¬ 
tion evaluations each iteration. The aim of this paper is to modify third order 
Newton-Secant method for solving nonlinear equations for multiple zeros with 
same order of convergence and efficiency index. In addition, the efficiency in¬ 
dex of an iterative method of order p requiring k function evaluations per 
iteration is defined by E(k,p) = ffp, see [5j. 

Let a is multi roots of fix) = 0 with multiplicity m and x a, then 
we can write f{x) = (x — a) m q(x ), where linu^a q(x) ^ 0. The function 
/ € C m [a , b] has a zero of multiplicity m at a in ]a, b[ if and only if 0 = /(a) = 
f'(a) = f"(a) = ■■■ = /( m_1 )(a), but ^ 0. Newton method for 

finding a simple zero a have been modified by Scheroder to find multiple zero 
of a non-linear equation which is of the form x n +i = x n — m converges 

quadratically m- 

Moreover, there have been many attempts to construct methods for finding 
multiple roots e.g. the work of Chun et al. [2], Dong [0], Hansen and Patrick 
[5], Osada (8], Victory and Neta JTj are proposed various iterative methods for 
finding multiple zeros a of a nonlinear equation /( x) = 0 where multiplicity 
to is known. 

This paper is organized as follows: Section [2] is devoted to the construc¬ 
tion and convergence analysis of a new method with convergence order three. 
Computational aspects, comparisons and dynamic behavior with other meth¬ 
ods are illustrated in Section [3] Finally, a conclusion is provided in Section 

H 


2 Description of the method and convergence analysis 


In this section, we propose a new modification of Newton-Secant’s method to 
find multiple zeros. Newton-Secant’s method is 


Vn — X n 
Xn +1 “ X n 


f(x n) 

/'(an) ’ 

_ /(a n) f(x n) 

f(.Xn)-f(.Vn) /'(an) ’ 


(n = 0,1,...), 


( 1 ) 


that order of convergence is three for simple roots. We aim at extending the 
method CD for multiple roots and build a method according to CD without 
any additional evaluations of the function or its derivatives by using of the 
parameter. In other words, convergence order of Newton-Secant method for 
approximating simple zero of nonlinear equations is three, while convergence 
order of this method is linear for finding multiple zeros. Therefore, we use a 
parameter 6 for solving this problem in the second term. 

We have 


Vn — %n 


/(an) 

/'(an) ’ 

_ Sf(x n ) f(Xn) 

9/(a»)-/(9n) /'(a n ) ' 


Xn+1 — X n 


( 2 ) 
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The order of convergence of the preceding method will be analyzed and the 
method will be adjusted accordingly to prove the following expected theorem. 


Theorem 1 Let a £ D be a multiple zero of a sufficiently differentiable func¬ 
tion f : D C R —> R for an open interval D with the multiplicity m, which 
includes Xq as an initial approximation of a. Then , the method fl|) has order 
three and 9 = (=l±m.y 1+m . 

Proof Let e n := x n - a, := y n - a, Ci := ^(1)^ ■ Using the fact 

that f(a) = 0, Taylor expansion of / at a yields 

f(xn) = e™ (c 0 + cie n + c 2 e 2 n + c 3 e 3 n ) + 0(4), (3) 

and 

/ (x n ) = e™ _1 (m + (m + l)cie„ + (m + 2)c 2 e 2 + (m + 3 )c 3 e 3 n + 0(4))- (4) 


Therefore 

f(xn) 1 ci 2 -(1 + m)c 2 + 2mc 0 c 2 3 , 4 . 

r/( \ 2 3 2 d" 

j (x n ) m m z co ui^Cq 

and hence 


(5) 


= y n -a = 


— 1 + m 


"On 


Cl 


IV? Cq 


— (1 + m)c 2 + 2mcoc 2 3 

- 


'0(4)- (6) 


For f{y n ) we also have 



fiVn) — (co + ci e,n,y + c 2 e 2 n . y + c 3 e 3 ny ) + 0(4,j,)- 

(7) 

Substituting @-0 in (fl5l). we obtain 



e n +i — D\e n + -D 2 e 2 + -D 3 4 + 0(4), 

(8) 

where 

e 

1_1+ 0+(- 1+m ) m rn’ 

1 \ m / 

(9) 

and 

9m~ 2+rn (—m{— 1 + m) m + 9m m (—1 + m))c 3 
(— 1 + m)((—1 + m) m — 9m m ) 2 co 

(10) 

and 

Om~ 3+m A 

(11) 


J 2(—1 + m) 2 ((—1 + m) m — 9m m ) 3 CQ ’ 


where 


A = (—1 + m) 2m (— 1 + m + 2 m 2 ){mc\ — 2(—1 + m)coC 2 ) 
+20 2 (—1 + m) 2 m 2m ((l + m)c 2 — 2mcoc 2 ) — 0(—1 + m) 1+m 
(m(3 + 4m) c? + 2(1 + m — 4m 2 )coC2). 


(12) 















4 


Massimiliano Ferrara et al. 


Therefor, to provide the three order of convergence, it is necessary to choose 
Di = 0 (i = 1, 2 ), so we have 


e = 


— 1 + to\ 1+ 
TO J 


(13) 


and the error equation becomes 

„2 


(i4> 

and method © has convergence order three, which proves the theorem. 


3 Numerical performance and dynamic behavior 

3.1 Numerical results 


In this section we apply the new method © to several benchmark examples 
and compare them with existing methods that have the same order of conver¬ 
gence. 

The new method is given by 


Un — %n 


f(x n ) 


__ ( — 1+rn) 1 + m f(Xn) f(x n ) 

x n+ i -x n - ( _ 1 +m )-i+m /(a . n) _ m -i+ m/(j/n) • 


The Osada’s method [ 8 ], is given by 


^n+1 — %n 


-to(to + 1) 


f{Xn) 

f(x n ) 



1 ) 


2 f'(Xn) 

f"{x n y 


(15) 


(16) 


The Dong’s method [4] , is given by 


Un = X n + 


f(x 


f'(x n 


Xn+l—Vn m( 1 J^) 1 m 

The Chun’s method [2], is given by 


„ _ _ _ m (( 2 7 ~l)"»+ 3 ~ 27 ) f(x„) 

Xn+1 — X n 2 /'( x „) 


7 (rra-l) 2 /'(an) 

2 f"(x n ) 


(1-7 )m 2 f(x n ) 2 f"{x n ) 

2 ■ 


(17) 


(18) 


In the numerical experiments of this paper we will use 7 = — 1. 

We have tested the method (ED on a number of nonlinear equations. To obtain 
a high accuracy and avoid the loss of significant digits, we employed multi- 
precision arithmetic with 100 significant decimal digits in the programming 
package of Mathematica 8 [B]. 

In order to test our proposed method m also compare it with the methods 
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ed, (ED and ED we compute the error, the computational order of conver¬ 
gence (COC) by the approximate formula (TTH 


coc 


In | (x n +i - a)/{x n - a)| 
In \(x n - a)/(x n -i - a)|' 


(19) 


And the approximated computational order of convergence, (ACOC) by the 
formula [3 


ln|(a;„ + i - x n )/{x n - x n -i) 


In | {x n 

- X n -i)/ (x n -i - X n - 2 ) | ' 

Test function f n 

root a 

fl(x) = (siiA x + x) b 

0 

h(x) = (In x + s/x - 5) 3 

8.3094326942315717953469556827 

f 3 (x) = (e x2 + 7x ~ 30 - l ) 6 

3 

f 4 (x) = (y/x- i - 3) 7 

9.6335955628326951924063127092 


Table 1 Test functions /i, ..., and root a. 



Method (1151) 

Method (1161) 

Method 1171) 

Method 1181) 

/l, x 0 = 0.1 

\xi - x*\ 

0.739e—3 

0.162e—2 

0.955e—3 

0.112e—2 

\X2 ~ X*\ 

0.364e—9 

0.106e—7 

O.llle—8 

0.215e—8 

X 3 - X* j 

0.434e—28 

0.304e—23 

0.175e—26 

0.149e—25 

COC 

3.0000 

3.0000 

3.0000 

3.0000 

ACOC 

2.9999 

2.9994 

2.9998 

2.9998 

f 2 , x 0 = 8 

\xi - x*\ 

0.166e—4 

0.719e—4 

0.277e—4 

0.376e—4 

\X 2 ~ X*\ 

0.246e—17 

0.876e—15 

0.189e~16 

0.667e—16 

x 3 - x*\ 

0.793e—39 

0.793e—39 

0.793e—39 

0.793e—39 

COC 

3.0000 

3.0000 

3.0000 

3.0000 

ACOC 

3.0000 

3.0000 

3.0000 

3.0000 

/ 3 , xo = 3.1 

\x\ - x*\ 

0.128e—1 

0.268e—1 

0.212e—1 

0.216e—1 

\X2 ~ X*\ 

0.325e—4 

0.105e—2 

0.333e—3 

0.362e—3 

x 3 - x*\ 

0.531e~12 

0.871e—7 

0.143e—86 

0.190e—8 

COC 

3.0000 

2.9964 

2.9993 

2.9993 

ACOC 

3.0004 

2.9428 

2.9857 

2.9858 

/ 4 , x 0 = 9 

\xi - I*| 

0.223e—3 

0.792e—4 

0.122e—4 

0.987e—4 

\X2 ~ S*| 

0.864e—14 

0.139e—15 

0.742e—19 

0.340e—15 

x 3 - x*\ 

0.918e—39 

0.918e—39 

0.918e—39 

0.918e—39 

COC 

3.0000 

3.0000 

3.0000 

3.0000 

ACOC 

3.0000 

3.0000 

3.0000 

3.0000 


Table 2 Errors, COC and ACOC for methods (115ft . (1161) . 1171) and (1181) . 
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In Table [2] our new method m is compared with the methods (UnD-dUD and 
(1181) on four nonlinear equations which were presented in Table [T| 


3.2 Dynamic behavior 

We already observed that all methods converge if the initial guess is chosen 
suitably. We now investigate the stability region. In other words, we numeri¬ 
cally approximate the domain of attraction of the zeros as a qualitative mea¬ 
sure of stability. To answer the important question on the dynamical behavior 
of the algorithms, we investigate the dynamics of the new method and com¬ 
pare them with common and well-performing methods from the literature. For 
more details one can consult mmmm- 

Let G : C — > C be a rational map on the complex plane. For z £ C, we 
define its orbit as the set orb(z) = {z , G{z ), G' 2 (z), ...}. A point Zq £ C is 
called periodic point with minimal period m if G m {z o) = zo , where m is the 
smallest integer with this property. A periodic point with minimal period 1 
is called fixed point. Moreover, a point zo is called attracting if |G'(zo)| < 1, 
repelling if |G'(zo)| > 1, and neutral otherwise. The Julia set of a nonlinear 
map G(z), denoted by J(G), is the closure of the set of its repelling periodic 
points. The complement of J(G) is the Fatou set F(G), where the basin of 
attraction of the different roots lie [T. 

For the dynamical point of view, in fact, we take a 512 x 512 grid of the 
square [—3, 3] x [—3, 3] £ C and assign a color to each point zq £ D according to 
the root to which the corresponding orbit of the iterative method starting from 
zq converges, and we mark the point as black if the orbit does not converge to 
a root, in the sense that after at most 100 iterations it has a distance to any of 
the roots, which is larger than 10~ 3 . In this way, we distinguish the attraction 
basins by their color for different methods. 

We have tested several different examples, and the results on the perfor¬ 
mance of the tested methods were similar. Therefore, we report the general 
observation here for following test problems which are presenter in Table [3] 


Test problem roots 

pi (A) = (z 3 - l) iu 1, - 0.5 ± 0.866025i 

p 2 (z) = (z 5 - z 2 4-1) 15 -0.808731, - 0.464912 ± 1.07147*, 0.869278 ± 0.388269* 

p 2 (z) = (2z 4 - z) 8 _0, - 0.39685 ± 0.687365*, 0.793701 


Table 3 Test problems pi,p 2 ,P 3 and their roots 
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Fig. 1 Comparison of basin of attraction of methods (U5i . net , inn and GU for test 
problem pi(z) = (z 3 — l) 10 



Fig. 2 Comparison of basin of attraction of methods CE3, GSJ. JED and GHJ for test 
problem P 2 {z) = (z 5 — z 2 + l) 15 



Fig. 3 Comparison of basin of attraction of methods inn , inn , inn and (1181) for test 
problem P 3 (z) = (2 z 4 — z) 8 


In Figures Ql [2] and El basins of attractions of methods (fl5l) . (fTfil) . (fTTl) and 
m with three test problems are illustrated from left to right respectively. As 
a result, in Figures [HE] basins of attraction of method (fl5l) seem to produce 
larger basins of attraction than other methods. 


4 Conclusion 

In this paper, Newton-Secant’s method for simple zeros was modified for find¬ 
ing multiple zeros of non-linear equations with same order of convergence and 
without any additional evaluations of the function or its derivatives. A numer¬ 
ical comparison with other methods shows that our new method is a valuable 
alternative to existing methods. In addition, a numerical investigation of the 
basins of attraction of the solutions illustrated that the stability region of our 
method it typically larger than that of other methods. 
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